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We present a high-precision numerical study of 3D random percolation viewed as a confining 
gauge theory. Using large correlation matrices among multiform Wilson loops we determine the 
low-lying masses in various spin channels. 
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1. The model: how to interpret pure percolation as a gauge theory 

In this Section we try to summarise the main points that make percolation theory a suitable 
framework in which to set up a gauge theory. The point of view from which is convenient to work 
is quite different from the standard one used in percolation theory (that mainly takes into account 
geometric aspects of the system), as pointed out in Ref. and is also well related to the two 
standard confinement mechanisms considered in gauge theories. 

The formulation of a gauge theory in terms of percolation is in some sense the most trivial 
gauge theory ever defined (its gauge group is the identity alone, G = {e}, thought for example as the 
^ — > 1 limit of a <7-state Potts model), but nevertheless contains all features of a satisfying confining 
gauge theory even if the model has no dynamics and every object in the system is completely 
independent. 

The ensemble of the model is described as follows: consider a lattice A (that in the following 
will be three-dimensional simple cubic), whose links are initially unoccupied. Then, every link can 
be occupied with a probability p and independently from all other links^. Thus, every configuration 
is some subset G of the lattice links. The quantity p plays the role of the coupling constant of the 
model, since it determines the mean size of the connected components {clusters) that constitute the 
configuration. In particular, there exists a well-defined value p = pc at which an infinite connected 
network of occupied links suddenly appears, called percolation threshold, and it is in fact a proper 
second order transition point with respect to the quantities that can be defined in the model, such 
as the cluster mean radius. 

The main observables in the theory are the Wilson loop associated to closed paths 7 on the 
dual lattice A: they are defined according to the rule 

Wy{G) = if G is hnked to 7 ; Wy{G) = 1 otherwise 

Note that such a definition does depend only on the subset BgCG of links belonging to loops, that 
is, observables do not change their value if dangling ends or "bridges" between loops are addded 
or removed. This can be seen as a sort of gauge-invariance of the theory, formulated in purely 
topological terms. 

In this formulation, naturally connected to ordinary gauge theories, the deconfined phase is 
identified with the disappearance of the infinite cluster, the transition point being mapped to a finite 
critical temperature. There are a number of numerical results that support the fact that the model, 
although having been drastically simplified, is a good description of confinement physics. 

2. Correlator functions and pure gauge mass spectrum 

It is typical for a pure gauge theory to possess a spectrum of physical states with well- 
defined quantum numbers, usually interpreted as bound states of gluonic degrees of freedom. The 
aim of this work is an accurate investigation on the glueball spectrum in the simple cubic three- 
dimensional bond percolation model. 

'We refer to that paper, and references therein, for a more detailed discussion on this introductory subject. 
^What is described here is the so-called bond percolation, but the situation is the same for the site percolation, apart 
from some implementation nuisances. 
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If one considers the correlator between plaquettes, 

'^(x,j) = (□,□,) -(□,>(□,) 

and then projects out nonzero momentum states, the result should encode all the sector of the mass 
spectrum carrying the plaquette quantum numbers = 0+ : 

(yi=t\)(xi=to) 

y X n=Q 

The lowest mass, which dominates for large t, is the inverse correlation length of the system. 
It has the same value found in the exponential decay of what is usually defined in percolation, the 
point to point correlator: 



G{x,y) = \ifx and y belong to the same cluster of ( 



G{x,y) = otherwise 



but the zero momentum projection of this object shows a single exponential law, and in addition 
the definition given here is not gauge-invariant, in the sense specified above. Thus, the former 
correlator couples with the whole spin/parity "channel" with vacuum quantum numbers, while the 
latter couples only to its lowest mass and is not suitable to analyse the entire spectrum, that could 
reasonably be constituted by an infinite tower of states. 

Due the the geometry of the model, the subgroup of rotations to consider is the dihedral group 
D4. As a consequence, there is no charge conjugation, and angular momentum is defined mod 4. 
Moreover, D4 has five irreducible representations such that the spin/parity families used for the 
classification of massive states are: 



0+ 


0- 


2+ 


2- 


1/3 



(notice the odd spin glueballs are grouped together and do not allow to determine parity). 

In an attempt to improve overlap with the actual glueball wave function, that is unknown, we 
extended the plaquette object to a basis of 0+ operators, from which we constructed the cross- 
correlation matrix: 



{y-xh=t 
I 



(o+)\ 



in which the subtracted term is nonzero only when working with the vacuum quantum numbers, 
i.e. in the 0+ channel. Our basis was made by the following 17 loops (called tetrises for short)^: 



?(o+) 



□ □ □ □ □ 



The same process can be applied to the other spin/parity channels; in these cases, however, one has 
to choose loop shapes not completely invariant under dihedral symmetry, and then construct proper 
linear combinations of the differently oriented copies of them. This approach has been successfully 
applied to the 3D Ising model spectrum in [0]. Example of observables constructed this way are: 



■'This definition includes also loops made by disconnected components, but treated as they were, i. e. there is no 
linking also if the cluster pierces (in opposite directions) and "ties together" two components. 
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Following this prescription, we constructed, for each non-0+ channel, a basis of operators by 
making use of the following tetrises: 



T- ^ 



*+*} 

}{. 11} 



For a given channel 7^, by diagonalising the cross-correlation matrices 'i^^ \t) one can iden- 
tify (effective) masses in the spin/parity family: this can be achieved with a naive diagonalisation 
for each value of f, or with a generalised eigenvalue problem (see, for instance, [01) by fixing a 
suitable in: 

^(f >fo)x = A'''(0^(fo)x 
Masses are then given, in inverse lattice spacings, by looking for a plateau in the limit 

Ut) 



rui = lim 



log 



A,-(f + l^ 



3. The Monte Carlo simulation: algorithm and setting 

We used lattices of size 60 x 60 x 100 (last direction was regarded as "time") with periodic 
boundary conditions. The critical percolation probability for such a lattice is p^^ ~ 0.248813, so 
we studied the confining range 0.256 < p < 0.262, where the correlation length does not exceed 
the value of about six lattice spacings. 

The algorithm is structured in such a way to examine a given configuration once for all mea- 
surements on a time-slice: first, a random configuration is constructed from scratch; then all dan- 
gling ends are removed from the configuration. From this reduced graph, for each time-slice t a 
table is constructed containing the values of all tetrises in all orientations and (summing over) all 
spatial positions; from these tables, operators' linear combinations are evaluated and eventually the 
crosscorrelation matrix. 

The construction of the (zero-momentum projected) table containing the value of each tetris 
in each orientation goes as follows. A first cluster reduction procedure is performed, but ignoring 
all time-like links passing through the (dual) surface t; in this way, the configuration is mapped to 
a list of associations between (occupied) links on i and two cluster labels, whose exact shape is 
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no more considered. Then, scanning links on t, and keeping track of the winding numbers while 
attaching clusters as prescribed by the mapping^, loop-like structures can be detected and their 
Unking with all tetrises can be checked at once. An alternative algorithm to perform the task of 
measuring topological linking of clusters has recently been proposed in [^]. 

On a modern one-CPU machine, processing a single configuration (which means considering 
more than a hundred different tetris shapes, times 60^ x 100 each) takes about two minutes of 
computation time. To have acceptable statistics, at least ~ 10^ configurations are needed. For each 
value of p we studied, the number of configurations generated was: 



0.256 


0.258 


0.260 


0.161 


52.000 


160.000 


52.000 


112.000 



4. Results: lightest scalar glueball 



From the highest eigenvalue in the 0+ channel, that is also the one showing the softest expo- 
nential decay, the mass of the lightest glueball is found. Its value is expected to scale as predicted 
by the 3D percolation critical index Vt,d — 0.8765 (whose value is known only numerically, see 



Oh- lightest mass 



m';{p)=Mt{p-p,) 



V3D 



p 


0.256 


0.258 


0.260 


0.262 




0.190(7) 


0.226(10) 


0.264(16) 


0.290(17) 



0.32 I 

0.3 - 
0.28 - 
0.26 - 
0.24 - 
0.22 - 
0.2 - 
0.18 



0.255 0.256 0.257 0.258 0.259 0.26 0.261 0.262 0.263 



Due to its low mass, the signal is well recognised up to ? ~ 15, while all other masses can 
be seen only in the first four or five lattice spacings. The results for wiq^ follow the expected be- 
haviour, and the scahng amplitude Mq^ = 13.31 ±0.43 is obtained, slightly larger than the estimate 

presented in [|I|]. The universal ratio ~ 4.46 is then evaluated. Its value is surprisingly close 
to the ~ 4.7 reported for SU{1) in the same dimensionality in [Q] and refined in [||], and of the 
same order of magnitude as the amplitude obtained for the Ising model. This fact confirms that the 
essential mechanism responsible for confinement is well included in the simpler percolation model. 

To give an estimate for this glueball's radius, square Wilson loops of side ranging from 1 to 15 
are evaluated, looking for the maximum coupling with the state: 



Squares' coupling to lowest state, p-0.258 



600000 — 

500000 - 

„ 400000 - 

g- 300000 - 

200000 - 

100000 - 

— 




((«) ~ 11 lattice spacings ~ 0.24 fm 
assuming physical string tension CT ~ (440MeV)^ 



6 8 10 
Square side 



To store the abstract graph, whose nodes are the found clusters, and its connections, a pointer based, tree-like 
standard structure is used. 
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5. Results: spectrum properties 



Following the same pattern, we obtained the lowest mass in each other channel, but somewhat 
less precisely because the noise drowned the signal already at f ~ 4-5 (the table below refers to 
p = 0.258): 



0+ 


0.226 ±0.010 





1.739 ±0.092 


2+ 


1.127±0.051 


2- 


1.131 ±0.079 


1/3 


2.417 ±0.378 



Lowest masses by chiinncl, p=0.258 



2.5 - 

2 - 



The spectrum presents striking resemblances to usual theories with nontrivial gauge symmetry: 
states in the families 2+ and 2^ well satisfy the expected degeneracy, and moreover, the proper 
hierarchy is found among channels: 



m, 



0+ 



< m,! < win < ni, 



1/3 



In these non-vacuum spin/parity channels, it is difficult to assign scaling amplitudes to lowest 
states: this is due to their apparent bad scaling behaviour, as well as their high mass damping the 

signal very soon. Furthermore, when looking for excited states in any channel (including 0+), the 

situation gets even worse since another issue comes up: mass estimates show a strong dependence 

on the choice of operators and their size. A detailed investigation on this problem, that could be a 

finite size issue, is currently under development. 

6. Conclusions 

This work has investigated on the percolation model a new aspect of what is expected from 
a reliable gauge theory, that is the presence of a physical spectrum of states with different mass 
and quantum numbers: despite some numerical difficulties encountered, the model is in agreement 
with what coud be expected. For more accurate results, a further analysis could be carried on, 
considering also the interest that this aspect of the percolation theory in itself could gain as well. 
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